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ABSTRACT 

In order to explain the slow rotation observed in a large fraction of accreting pre-main-sequence 
stars (CTTSs), we explore the role of stellar winds in torquing down the stars. For this mechanism 
to be effective, the stellar winds need to have relatively high outflow rates, and thus would likely be 
powered by the accretion process itself. Here, we use numerical magnetohydrodynamical simulations 
to compute detailed 2-dimensional (axisymmctric) stellar wind solutions, in order to determine the 
spin down torque on the star. We discuss wind driving mechanisms and then adopt a Parker-like 
(thermal pressure driven) wind, modified by rotation, magnetic fields, and enhanced mass loss rate 
(relative to the sun). We explore a range of parameters relevant for CTTSs, including variations 
in the stellar mass, radius, spin rate, surface magnetic field strength, the mass loss rate, and wind 
acceleration rate. We also consider both dipolc and quadrupole magnetic field geometries. 

Our simulations indicate that the stellar wind torque is of sufficient magnitude to be important for 
spinning down a "typical" CTTS, for a mass loss rate of ~ 10~^Mq yr~^. The winds are wide-angle, 
self-collimated flows, as expected of magnetic rotator winds with moderately fast rotation. The cases 
with quadrupolar field produce a much weaker torque than for a dipole with the same surface field 
strength, demonstrating that magnetic geometry plays a fundamental role in determining the torque. 
Cases with varying wind acceleration rate show much smaller variations in the torque suggesting that 
the details of the wind driving are less important. We use our computed results to fit a semi-analytic 
formula for the effective Alfven radius in the wind, as well as the torque. This allows for considerable 
predictive power, and is an improvement over existing approximations. 

Subject headings: accretion, accretion disks — MHD — stars: magnetic fields — stars: pre-main- 
sequence — stars: rotation — stars: winds, outflows 



1. INTRODUCTION 

For more than half a century, the spin rates and the 
angular momentum evolution of stars have been topics 
of vigorous study. We know that stellar winds are re- 
sponsible for the spinning down of late-type (later than 
F2) main sequence stars (Parker 1958; Schatzman 1963 
Kraft 1967; Skumanich 1972; Sodcrblom 1983; Kawaler 
1988i; iMacGreeor fc Brenner 1991, ; Barnes fc Sofia, 1996 ; 
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jjQuvier et al.llltjijYI ). There is still progress t o be made 
on main sequence star spins (|Barnesl l2003h . but per- 
haps the largest open questions remain at the pre-main- 
sequence phase, which determines the "initial condi- 
tions" for the spin histories of stars. 

By the time intermediate/low mass (< 2Mq) pre- 
main-sequence stars become optically visible (T Tauri 
stars; TTSs), they already have ages around 10^ - 10^ 
yrs. A large fraction of TTSs (called classical TTSs; 
CTTSs) are observed to actively accrete material from 
a dis k at a rate within a wide ran ge of ~ 1Q~^Mq yr^^ 
(e.g.. lJohns-Krull fc Gaffordll2002l ). At this rate, the an- 
gular momentum accreted from the orbiting disk should 
spin up the stars to a substantial fraction of breakup 
speed in a short amount of time (comparable to their 
ages) . The fact that th e stars are also still contracting 
(e.g.. lRebull et al.ll2002f ). and that they presumably were 
accreting at much higher rates before they became opti- 
cally visible, further adds to expectation of fast rotation. 

Large data sets for the spins of TTSs in s tar formation 
regions and clusters of different ages (see iRebull et al.l 



l2004l for a compilation) show that approximately half 
of the stars are rotating rapidly and do seem to spin up 
as ex p ected a s they app roach zero - age main sequence 
JVogel & Kuhi 1981; B ouvier et all 119971 ; iRebull et all 
[2004; ,IIerbst et al.. 2007l). However, the surprise is that 
the other ~half of TTSs exhibit much slower rotation 
rates (~ 10% of breakup speed) at all ages. Recent 
studies have shown a correlation between slow rotation 
and the presence of a n a accretion disk (see especially 
ICieza fc Balibeij[2007() , though this ide a has b een contro- 
versial in the past (e . g..lStassun et an[l 999; H erbst et al.l 
120001 IStassun et all |2001| ; iHerbst et al.l l2002i ). This is 
still an open issue, but it is clear than an efficient angu- 
lar momentum loss or regulation mechanism is operating 
for the slow rotators. 

Although alternative ideas ha ve been proposed since 
(|Konig]||1991l ; iShu et al.lfT99l see iMatt fc Pudritgi2007bl 
for a historv). iHartmann fc Stau ffcr (1983) offered the 
first potential explanation for the slow rotators, namely 
that massive stellar winds may be responsible for 
carrying off substan tial angular momentum (se e also 
iTout fc Prin"ii3ll992f ). In iMatt fc Pudrit j (|2005al here- 
after Paper I), we extended this idea to consider the ef- 
fects of the magnetic interaction between the star and 
disk, and we used a 1-dimensional scaling from the so- 
lar wind angular momentum loss to estimate the torque 
for TTSs. The scaling suggested that, for an obser- 
vationally c onstrained dipole magnetic field strength of 
200 G (;e.g.. lJohns-Krull et al.lfl999l ; iBouvier et al.ll2007t 
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Johns-Krulll l2007allbl : ISmirnov et all l2003at lYang et all 
2007( ). it might indeed be possible for the stellar wind to 



extract enough angular momentum to explain the slow 
rotators. For stellar winds to balance the accreted an- 
gular momentum, the wind outflow rate needs to be a 
substantial fraction of the accretion rate. In Paper I, we 
suggested that this is possible, if a fraction of the energy 
liberated by the accretion process actually powers the 
stellar wind. 

The pre-main-scqucn ce phase is, in fac t, marked by 
powerful outflows (Reipurth fc Balivll2001h . In the most 
powerful sources, due to the large linear momenta of the 
outflows ("Konigl fc Pudrit j l2000l ). the X-ray luminosi- 
ties JP ccampli 198ll), and possible de tection of rotation 
jBacciotti ct al. 2002; Anderson et al. 2003: C offey et al.l 
[2004: Fcrreira_ct al. 2006; Coffcy_ctjil. 2007), it appears 
that most of the flow arises from the accretion disk, 
rather than the star. It is not clear what fraction of 
the total outflow may actually originate from the star, 
and thus how powerful are the stellar winds compared to 
main-sequence phase winds or compared to their accre- 
tion rates. 

There is some observational evidence for powerful stel- 
lar winds from CTTSs, as distinguished from inner disk 
winds. In particular, lEdwards et"al] (|2003l |2006[ ) ob- 
served the He 1 10830 A line in 39 CTTSs and saw several 
cases with a broad, deep, blue-shifted absorption, indi- 
cating outflow velocities of typically a few hundred, and 
up to ^ 400 km s^^. They concluded that this feature 
is best interpr eted as arising in a n optically thick stellar 
wind (see also lDupree et al.(l2005i) . They also suggested 
the winds may be accretion-powered, since the wind sig- 
natures are most prevalent in the stars with highest ac- 
cretion rates and absent in non-accreting sy stems. Sub- 
sequent modeling of the He 1 10830 A line bv' Kwan et al.l 
(|2007, ) indicates that approximately half of these CTTSs 
show evidence for a po werful stellar wind. Furthermore, 
[Kurosawa et al.l ()2006[ ) modeled the Ha emission line in 
these systems and suggested that a stellar wind compo- 
nent could most naturally explain the profiles obse rved in 
^ 7% of the stars in a sample compiled by Reipur th et al.l 
([1991). 

There already exists some theoretical work on stel- 
lar winds, specifically from pre-main-seque nce stars, 
with a focus o n the wind dri v ing m echanism (|Decamplil 
119811: ,Hartm ann~et al.l Il982l Il990l) or the collimation 
of th e winds~ ([Fendt et aiTliggSl : iFendt fc Camenzindl 
|1996[) . These do not discuss the expected angu- 
lar momentum outfiow rates, however. The works 
that do calcula te stellar wind torques for pre-main- 
sequence stars (|Hartmann fc Ma cGrcgor 1982: Mestelj 
1981 iHartmann fc Stauffcr' 119891 : iTout fc'Pringltj figgT 



Paatz fc Camenzind 4996. , Paper I) are either based on 



a 1-dimensional formulation and/or have made a priori 
simplifying assumptions regarding the stellar magnetic 
field structure, wind fiow speed, and latitudinal depen- 
dence of the wind. Calculating the stellar wind torque 
reliably is a complex, multi-dimensional problem, and 
more work is needed to develop the stellar wind theory 
further. 

The primary goal of this paper therefore, is to take 
the next major step in developing the accretion-powered 
stellar wind picture by rigourously computing the steady- 



state solutions of winds from spinning magnetized stars. 
We carry out a parameter study to provide a range 
of possible solutions that are expected to charaterize 
accretion-powered stellar winds. Where possible, we 
compare our results to analytic magnetohydrodynamic 
(MHD) stellar wind theory. In a companion paper, 
we will use these solutions to compare the stellar wind 
torques and wind driving power with the torque and en- 
ergy deposition expected to arise from the interaction of 
the star with its accretion disk. 

In the following section (^24}, we give a brief intro- 
duction to basic stellar wind theory. This provides the 
motivation for using a numerical approach and sets the 
stage for comparing our numerical results with the an- 
alytic theory. Section 12.21 contains a discussion of our 
adopted wind driving mechanism. We describe our nu- 
merical method for obtaining solutions in section [31 and 
present the results in section[4l Section[5lcontains a semi- 
analytic formulation for the torque and a comparison to 
previous theory. 

2. MAGNETIZED STELLAR WINDS: NEEDED 
BACKGROUND 

2.1. Magnetic Stellar Wind Theory 

Standard MH D wind theory (i .e ., mag netic rotator the- 
ory), following [Weber fc David (|1967[ ). characterizes a 
steady-state flow of plasma along a magnetic field line 
that is anchored to a rotating object, which we will here- 
after take to be a star. One of the key results is that 
the angular momentum outfiow r ate per unit mass loss 
is given very simply as (s ee, e.g., [Weber fc Davia4967[ : 
IMestell[l96a IMichellll969l ) 



I — rj^r^. 



(1) 



where fi* is the angular rotation rate of the star, and 
TA is the cylindrical radius at which the outfiow speed 
equals the local magnetic Alfven speed. 
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where p is the local mass density and Bp is the strength 
of the poloidal magnetic field. Bp ~ {B^ + i?^)^/^, in 
cylindrical (r, cj), z) coordinates. Equation ([1]) indicates 
that the quantity of angular momentum carried in the 
wind is as if the wind material is corotating out to r^ 
and conserves its angular momentum thereafter. Thus, 
TA is often referred to the magnetic "lever arm." In re- 
ality, the azimuthal velocity of the wind, v^, is a smooth 
(i.e., differentiable) function of radius, and the difference 
between Vipr and I at all radii equals the torque trans- 
mitted by azimuthally twisted magnetic field lines. 

By integrating the mass flux times I over any surface 
enclosing the star, one obtains an expression for the total 
angular momentum outflow rate and, by Newton's third 
law, the torque on the star: 



(3) 



where is the integrated wind mass loss rate. Since 
the value of ta will generically not be the same along 
each field line, equation ([3|) defines the quantity (rA), 
wh ich is the mass-loss- weight ed average of r\ (suggested 
bv [Washimi fc Shibatal[1993l ). Hereafter, we will simply 

1 /2 

refer to this average as ta = (r^) 
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The difficulty now lies in calculating ta- The lever 
arm length clearly depends on the stellar surface field 
strength (B*), stellar radius (i?*), and Mw because these 
directly affect Alfven condition. But it also depends on 
the flow speed and field structure, which are not possi- 
ble to determine a priori in the wind. The flow speed is 
influenced by the thermal energy in the wind as well as 
rotation. In ad dition, there exist two different regimes 
(jBelcher &: Ma cGregor 1976): the fast magnetic rotator 
regime, where the flow speed is mostly determined by 
magnetorotational effects; and the slow magnetic rotator, 
where the flow speed is solely determined by the wind 
driving. The field structure in the wind, even though 
the geometry may be known at the stellar surface, is de- 
termined by the self-consistent interaction between the 
wind and rotating magnetic field and thus is a function 
of all parameters. Therefore one can only calculate ta 
by making a priori assurngtions ab out the field struc- 
ture and/or flow speed jW eber & David Il967t iMestel 
1968': OkamotdllQT-t [Mestel 1984; M estel fc SDruitlll987l : 
Kawaler 1988 ) or by using iterative techniques (or numer- 
ical simulations; iPneuman fc Kor)dll971t lSakurail[l985|: 



Washimi fc Shibata| 119931 : iKeppens k Goedbloedl boOO : 
Matt & Balick 200|). 

All of these methods are complementary. The analyt- 
ical work, in which the held structure is guessed, pro- 
duces a predictive f ormulation of the stellar wind torque 
(e.g.. lKawaleijl988l ). However, usually the formulation of 
the field structure introduces more parameters (such as 
a power law index for the magnetic field) , so that almost 
any result can be obtained by adjusting these. Further- 
more, the field structure in the analytic models has no 
explicit dependence on (e. g.) r2^,, whic h is exhibited in nu- 
merical simulations (e.g. jMatt fc Ba lick 2004). The nu- 
merical simulation technique has the advantage of calcu- 
lating the field structure and flow speed self-consistently. 
However, a single simulation does not predict the de- 
pendence of TA on parameters, and to date, not enough 
parameter space has been explored. Thus, to date, there 
exists no formulation for the stellar wind torque that con- 
vincingly applies over a wide range of conditions (e.g., 
over a range of i?,, Af„, and fi*). 

In this paper we will use 2-dimensional (axisymmet- 
ric) MHD simulations to calculate the torque and cor- 
responding value of ta. This will allow us to check the 
estimate for ta of Paper I (and previous works) . In ad- 
dition, we will carry out a parameter study to determine 
the dependence of the stellar wind torque on parame- 
ters, over a range of conditions appropriate for TTSs, 
and compare with the predictions of analytic theory. 

2.2. Wind Driving Mechanism 

It is not known what dr ives winds from TTSs. These 
stars have act i ve co ronae (Feigelson fc Montmerldll999l : 
IStassun et all 120041 : Havata et al. 2005), and it thus 
seems a reasonable assumption that they also drive solar- 
like coronal winds in which thermal pressure plays a sig- 
niflcant r ole in the wind acceleration. Based o n a calcula- 
tion fr om iBisnovatvi-Kogan fc LamzinI (|1977t ) , iDecamplil 
()1981h concluded that, in order for the wind emission to 
be consistent with the X-ray observations, the mass loss 
rate of a T Tauri star coronal wi nd must be le s s than 
~ 1O"^M0 yr^^ Furthermore, iDupree et all ()2005f ) 



found evidence for a stellar wind with a coronal tem- 
perature in the CTTS TW Hya (though this conclu sion 
has been challenged bv lJohns-Krull fc Herczedl2007l ). 

The assumption of thermal pressure driving is a sim- 
plification, even for the solar wind. It is known that a 
major factor in driving the solar wind is Alfven wave 
momentum and energy deposition. Two important re- 
cent studies have done self-consistent analyses of the 
combine d problem of both solar wind heating and acce l- 
eration (jSuzuki fc Inutsukall2006l : ICranmer et"alll2007f) . 
The first paper shows that low frequency, transverse mo- 
tions of open field lines at the photosphere leads to tran- 
sonic solar winds for superradial expansion of the wind 
cross-section. If the amplitude of these transverse pho- 
tospheric motions exceeds 0.7 km s~^, fast winds are 
produced and the dissipation of wave energy heats the 
atmosphere to a million degrees. The results are sensi- 
tive to the amplitude of the velocity perturbations, and 
the simulations show that the solar wind vitually dis- 
appears for amplitudes < 0.3 km s~^. These numerical 
simulations also show that Alfven wave pressure dom- 
inates the gas pressure in the solar acceleration region 
(1.5i?0 < R < IORq). The second paper shows sim- 
ilar results. This work shows that there are three key 
parameters that control wind heating and acceleration: 
the fiux of acoustic power injected at the photosphere, 
the Alfven wave amplitude there, and the Alfven wave 
correlation length (characterizing wave damping through 
turbulence) at the photosphere. 

Our primary goal here is to evaluate the angular mo- 
mentum transported away from the star by the stellar 
wind. Thus, in this work, we do not discuss the thermo- 
dynamic properties of the wind and instead focus on the 
angular momentum transport. Fortunately, this torque 
does not much depend on what drives the wind. Rather, 
the torque depends primarily on the stellar magnetic 
field, rotation rate, radius, M^, and the wind velocity. 
As long as "something" accelerates the wind to speeds 
similar to what we see in our simulations, the torque we 
calculate will be approximately correct. 

We expect that the Alfven waves in accreting TTS 
winds will have a significant, if not dominant contribu- 
tion to both the acceleration and heating of their winds. 
These waves will be launched along the open field lines 
that originate from the TTS photosphere at latitudes 
comparable to those that harbour field lines carrying the 
accretion flow onto the star. The irregular accretion flow 
should generate very large (i.e., much larger than acous- 
tic motions in the solar photosphere) acoustical trans- 
verse motions in the TTS photosphere as it impinges 
upon the star. These large amplitude perturbations, gen- 
erated by the accretion flow itself, may be the ultimate 
driver for the Alfven wave flux that drives our proposed 
accretion-powered stellar wind. 

Note that the driving force can be parameterized as be- 
ing prop ortional to — (where ^ is the wave energy den- 
sity; [Decanicli||l981l)- This has the same functional form 
as the thermal pressure force (—VP) used in our sim- 
ulations. Severa l authors (e.g.. iHartmann fc MacGregoH 
[19801: IDecamplil [19811 : iHolzer et all 119831 : iSuzukil l2007f) 
computed velocity profiles for cool (~ 10^ K) Alfven 
wave-driven winds. These works exhibit wind velocity 
profiles that are similar to what is expected from ther- 
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mal pressure driving of hotter winds. Therefore, we can 
think of thermal pressure driving as a proxy for some 
other driving mechanism. Also, it will be important to 
have these solutions to compare with future work that 
includes different driving mechanisms. 

In this paper, we restrict ourselves to mass loss rates 
of Mw < 2 X 1 0~^M(7) yr~^ As justified above, we adopt 
a Parker- like ( Parkeiiri958l ) coronal wind driving mech- 
anism, modified by magnetic fields, stellar rotation, and 
an enhanced mass loss rate (relative to the sun). As 
the nature (e.g., temperature) of TTS stellar winds is 
not well-known, our detailed solutions of coronal winds 
will enable us to look at the expected radiative prop- 
erties, a posteriori, allowing for further constraints on 
rea l systems. We will show in a forthcoming paper (and 
see iMatt fc Pudrit j l2007af) that the expected emission 
from the simulated winds presented here rules out ther- 
mal pressure driving at a subst antially lower mass loss 
rate than the limit of iDecamplil (1981). 

3. NUMERICAL SIMULATION METHOD 

We calculate solutions of steady-state winds from iso- 
lated stars (no accretion disk) , using the finite-difference 
MHD code of Matt & Balick (2004), and the reader wiU 
find further details there^ (and references therein). As- 
suming axisymmetry and using a cylindrical (r, 0, z) 
coordinate system, the code employs a two- step Lax- 
Wendroff scheme (jRichtmver fc MortonI Il967l ) to solve 
the following time-dependent, ideal MHD equations: 



dp 
dt 
d{pv) 
dt 
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where p is the density, v the velocity, P the gas pressure, 
G Newton's gravitational constant, M, the stellar mass, 
R the spherical distance from the center of the star {R^ — 
-\- z^), e the internal energy density, B the magnetic 
field, J the volume current, E the electric field, c the 
speed of light, and 7 the ratio of specific heats. 

To obtai n steady-state w i nd so lutions, we follow the 
method of iMatt fc Bahckl (l2004f). which is also simi- 
lar to that employed b v lWashimi fc Shibatal |1993.) and 
iKeppens fc Goedbloedl (|1999( ). It involves initializing 



1 IMatt fc BalTcH 1I2OOI) ran cases with isotropic hydrodynamic 
variables at the base of the wind and also cases with enhanced 
polar winds. Here we only consider the isotropic case. 



the computational grid with a sp herically sym metric, 
isothermal Parker wind solution (jParked Il958l ). plus 
force-free dipole (and sometimes quadrupole) magnetic 
field. When the simulation begins, the wind solution 
changes from the initial state due to the presence of the 
magnetic field, the rotation of the star, and the poly- 
tropic equation of state (P cx p^). The simulations 
run until the system relaxes into a steady-state (within 
a small tolerance) MHD wind solution. The code uses 
nested computational grids so that the wind can be eas- 
ily followed to large distances (several tens to hundreds 
of i?*). 

This method results in a steady-state solution for the 
wind that is determined solely by the boundary condi- 
tions held fixed at the base of the stellar corona (the 
"stellar surface" ) . In order to capture the appropri- 
ate physics within the framework of a finite difference 
scheme, we employ a four-layer boundary for the star, on 
which the various physical quantities are set as follows. 
We consider the spherical location i? = 30, in units of the 
grid spacing, to be the surface of the star. For all grid- 
points such that R < 34.5, the poloidal velocity is forced 
to be parallel with the poloidal magnetic field {vp \\ Bp, 
where the poloidal component is defined as the vector 
component in the r-z plane). Where R < 33.5, p and 
P are held constant (in time) at their initial values. For 
R < 32.5, tip is held at zero, while is held at coro- 
tation with the star. For R < 31.5, -Bp field is held at 
its initial, dipolar value, while is set so that there is 
no poloidal electric current at that layer (which gives it 
a dependence on the conditions in the next outer layer, 
31.5 <R< 32.5). 

These boundary conditions properly capture the be- 
havior of a wind accelerated from the surface of a rotat- 
ing magnetized star, as follows. There is a layer on the 
stellar boundary {R > 32.5) outside of which the velocity 
not fixed, but is allowed to vary in time. In this way, the 
wind speed and direction is not specified, but is deter- 
mined by the code in response to all of the forces. By 
holding P fixed at its initial value for all R < 33.5, we 
constrain the pressure gradient force (thermal driving) at 
the base of the wind to be constant in time. Also, hold- 
ing the density fixed at i? < 33.5 allows the region from 
where the wind flows to be instantly replenished with 
plasma. Thus, the base of the wind maintains a constant 
temperature and density, regardless of how fast or slow 
the wind flows away from that region. The existence of 
a layer in which Vp = and Bp can evolve (namely, at 
31.5 < R < 32.5) allows Bp (and Vp) to reach a value that 
is self-consistently determined by the balance of magnetic 
and inertial forces. We set the poloidal velocity parallel 
to the poloidal magnetic field for the next two outer lay- 
ers, to ensure a smooth transition from the region of pure 
dipole field and zero velocity to a that with a perturbed 
field and outfiow. Setting B^ so that the poloidal elec- 
tric current is zero inside some radius ensures that the 
field behaves as if anchored in a rotating conductor (the 
surface of the star). Also, this ensures that B^ evolves 
appropriately outside the anchored layer according to the 
interaction with the wind plasma. 

The key physical parameters can be represented by 
the characteristic speeds of the input physics, namely 
the sound speed at the base of the corona, Cg, the escape 
speed from the surface of the star, Vcsc, the rotation speed 
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of the star, and the Alfven speed at the base of the wind. 
We specify the ratio of Cs/wcsc as our parameter, rather 
than the sound speed alone. This seems the most reason- 
able, since the temperature of a thermally driven wind is 
regulated somewhat by the interplay between the ther- 
mal energy input and the expansion of the corona (the 
wind) against gravity. To first order, a hotter wind ex- 
pands more rapidly against gravity allowing less time for 
the gas to heat, and a cooler wind expands more slowly, 
allowing more time to heat. Once the value of the stellar 
mass and radius is specified, the ratio of Cg/wcsc deter- 
mines the temperature held fixed on the stellar boundary, 
as described above. The wind plasma is characterized by 
a polytropic equation of state, and so 7 is also a param- 
eter. We parameterize the stellar rotation rate as the 
fraction of breakup speed, 

/ = n,0\GAQ-'^\ (11) 

The Alfven speed is determined by the magnetic field 
strength and coronal density. Rather than taking the 
Alfven speed as a key parameter, we specify the field 
strength at the equator of the star {B^) as our parame- 
ter, in order to connect the simulations as much as pos- 
sible to observationally constrained quantities. For the 
same reason, we specify Mw as a parameter, rather than 
the coronal density. In the simulations, we must spec- 
ify the base density, p*, to be held fixed on the stellar 
boundary, and the value of in the steady-state wind 
is not solely determined by p*. For example, the rota- 
tion of the star can enhance via magneto-centrifugal 
flinging, and a strong magnetic field can decrease 
by inhibiting fiow from a region near the equator that 
remains magnetically closed (the "dead zone" ) . In other 
words, Mw is not an a priori tunable parameter; rather, 
it is a result of the simulations. Therefore, to treat Mw as 
our tunable parameter, we adopt an iterative approach. 
This entails first running a given simulation with a guess 
for p*, checking the resulting value of Afw, and then ad- 
justing p* and rerunning the simulation. We iterate until 
the desired value of M^ is achieved (within a tolerance 
of 2%). This typically required 2 to 4 iterations, so the 
ability to treat M^ as a chosen parameter comes at a 
substantial cost. 

4. STELLAR WIND SOLUTIONS 

4.1. The Fiducial Case 

We start by presenting the results of our stellar wind 
simulation for parameters with values that represent a 
"typical" T Tauri star and follow th e fiducial values of 
Paper I and lMatt fc Pudritd (|2005bf ). Table [T] lists the 
fiducial parameters. We consider a low mass pre-main- 
sequence star, with a surface escape speed of Vesc ~ 309 
km s~^. A dipole magnetic field stre ngth of 200 Gauss 
is co nsistent with 3g upper lim its (jJohns-Krull et aTl 

fig"; Smirnov ct al. 2004, 2003b) or marginal detection 
Qirnoy_et_ al. 2003a; Yang etal. 2007.) of the longitudi- 
nal magnetic field measured for CTTSs. We seek primar- 
ily to understand the slow rotators, for which a rotation 
rate of 10% of breakup is appropriate. In Paper I, we 
estimated that an accretion-powered stellar wind for a T 
Tau star mig ht have « 1.9 x IO'^Mq yr"!, so we 
use this as our fiducial value. 



TABLE 1 
Fiducial Stellar Wind 
Parameters 



Parameter Value 



M, 0.5 Mq 
2.0 R(7) 

B* (dipole) 200 G 

/ 0.1 

Mw"" 1.9 X IO-^Mq yr-l 

Cs/«osc 0.222 

7 1.05 



^ In order to treat Mw as a param- 
eter in the simulations, our method 
is to adjust the mass density at the 
base of the wind until the desired 
Mw is achieved in the steady-state. 
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Fig. 1. — Fiducial case: greyscale of log density, velocity vectors, 
and magnetic field lines illustrate the structure of the steady-state 
wind solution (see table [TJ. The dashed line represents the Alfven 
surface, where the wind speed equals the local Alfven speed. The 
rotation axis is vertical, and the longest vector corresponds to 160 
km . Black corresponds to a density above 5.3 X 10"^"^ g cm~^ 
and white to a density below 2.6 X 10~^^ g cm~^. 

In a thermally-driven wind, the coronal sound speed 
should be comparable to the escape speed, and we use 
Ca/vcsc = 0.222 as our fiducial value. This value gives 
wind speeds that are appropriate in the solar case. The 
choice of polytropic index 7 is also important. At large 
distances (~AU) from the sun, the solar wind plasma 
is well characterized by an effective 7 be t ween approxi- 
mate ly 1.5 and 5/3 (jFeldman et al.lfl998l : iKrasnopolskvl 
l2000f ). However, in the region where the wind is ac- 
celerated (within a few solar radii), thermal conduction 
and other heating a nd cooling effects play a role (e.g., 
ICranmer et al.l[2007f ). resulting in an effective 7 closer to 
unity (is othermal). Our fiducial va lue o f 7 = 1.05 was 
used by IWashimi fc Shibatal (|1993f ) and iMatt fc Bailed 
(|2004D for solar-like winds. This nearly isothermal value 
of 7 approximates the thermodynamics of a gas with a 
true value of 7 = 5/3 that is heated as it expands. 
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TABLE 2 

Stellar, Wind Torques and Lever Arm Lengths 



Case 


p. 


M„ 


Tw 


/ 2 \l/2 




(10-"^) 


, 9M0, 


(10^"^ erg) 


(R,) 


nducial 


3.67 


1.89 


1.77 


6.97 


/ = 0.004 


7.62 


1.86 




8.33 


/ = 0.2 


1.36 


1.87 


2.82 


6.26 


/ = 0.05 


6.01 


1.88 


1.06 


7.65 


£>t = 400 G 


3.67 


1.86 


3.27 


9.55 


B* = 2 kG 


o.u / 


1 92 


-LO.O 


1 Q 


1 kG quad. 


2.92 


1.87 


1.37 


6.17 


2 kG quad. 


4.38 


1.93 


2.11 


7.53 


low Mw 


0.377 


0.187 


0.500 


11.8 


very low Mw 


0.0755 


0.0378 


0.204 


16.7 


R* = 1.5 Rq 


5.71 


1.86 


1.10 


5.96 


R* = 3 R(7) 


1.99 


1.89 


3.43 


8.75 


M« = 0.25 M0 


5.06 


1.91 


1.47 


7.52 


M, = 1 Mq 


2.59 


1.88 


2.11 


6.42 


Cs/Vesc = 0.245 


0.773 


1.87 


1.59 


6.64 


Cs/'Uesc = 0.192 


55.4 


1.89 


1.91 


7.23 


7 = 1.10 


11.1 


1.87 


2.19 


7.79 
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Figure [T] shows the result of our fiducial case sim- 
ulation, which illustrates the steady-state wind solu- 
tion. The non-spherical shape of the Alfven sur- 
face (which eventually crosses the rotation axis at 
larger radii than shown) is m ainly due to magnetorota- 
tional effects in the wind (see lWashimi fc ShibatalllQQSl : 
iMatt h Balickll2004f ). This demonstrates that the fidu- 
cial T Tauri star wind exists in the thermo-centrifugal 
regime where thermal and magnetocentrifugal ef fects are 
of similar importance for accele rating the wind (jSakurail 
fl985l: iWashimi fc Shibatal fl993l) . These winds are self- 
collimated, while still exhibiting substantial fiow at all 
latitudes. 

From the simulation, we calculate Mw and the total 
angular mome n tum o utflow rate, t„, as described by 
IMatt fc Balicg (|2004D . Then, using equation ([SD, we 

1 /2 

calculate the effective lever arm length, ta = (/"a) 
These results are listed in the first row of table [2l where 
we also list the coronal base density p* that we iteratively 
chose to give the desired value of Mw 

The wind base density of ~ 10~^^ gm cm^'^ is 5 or- 
ders of magnitude lar ger than required for si mple solar 
wind models (e.g., iWashimi fc Shibatal Il993f) . This is 
expected, since the fiducial Mw is 5 orders of magnitude 
higher than the solar value, and the wind speeds are com- 
parable. 

The fiducial stellar wind torque of « 1.8 x 10^^ ergs is 
capable of balancing the spin up torque from accretion at 
a rate of 4.4 x 1O~^M0 yr~^. The basic conclusion here is 
that the stellar wind torque for the fiducial case is of the 
right magnitude to be important for spinning down the 
star, as required by the accretion-powered stellar wind 
scenario. We chose our fiducial parameters to compare 
with the estimate of Paper I that^ rp^jR^ k, 12.2. We 
can see that their estimate, based on scaling of ID wind 
theory from solar values, was a 75% overestimate of j'a- 
We will identify the reasons for this in section 14.2.11 

^ Paper I actually quotes a value of rj^/Rt, = 15, but our defi- 
nition differs slightly here (compare eq. [3]here with eq. 2 of Paper 
I), so the lever arm length corresponds to 12.2 i?« here. 



Fig. 2. — Same as figure [T] but for the / = 0.004 (similar to 
solar) case. The maximum velocity vector corresponds to 95 km 
s-i. 

4.2. Parameter Study 

To establish the dependence of ta on parameters and 
to calculate wind solutions that are applicable to a wide 
range of conditions that are observed or often assumed 
for T Tauri stars, we carried out a limited parameter 
study with our simulations. The results are listed in 
table [21 and we briefly discuss each case below. The first 
column in the table lists the value of the parameter that 
is changed relative to the fiducial case. For each case, all 
other parameters are identical to the fiducial case. Note 
that since we consider Mw as a key parameter, the value 
of varies from case to case. 

4.2.1. Spin Rate 

As with the fiducial case, all but one simulation in our 
parameter study lie in a regime that is near the bound- 
ary between slow and fast magnetic rotators. The one 
exception is a case with a fractional rotation rate equal 
to the solar value oi f — 0.004, which represents a slow 
magnetic rotator. Figurc[2]illustrates the structure of the 
steady-state wind solution for this case in the same for- 
mat as the previous figure. A comparison between the 
two figures reveals that rotation indeed infiuences the 
detailed structure of the velocity field and the magnetic 
field in the wind, which manifests itself as a difference 
in the shape of the Alfven surface. All else being equal, 
the effect of faster rotation is to reduce the effective lever 
arm length, as evident in table [2] (see also ISakurai|[l985t 
IWashimi fc Shibatal [1993( 1. Ahhough the qualitative ef- 
fect of rotation on the shape of the A lfven surface was an- 
ticipa ted in analytic theory (e.g.. IBelcher fc MacGrego^ 
Il976( l. this effect is not properly included in any existing 
analytic formulation for calculating the torque. This is a 
primary reason that numerical simulations are required 
to convincingly calculate the self-consistent wind solu- 
tion, especially when considering winds that exist near 
the boundary between slow and fast magnetic rotators. 
The effect of rotation on ta was not considered in the 
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Fig. 3. — Same as figure [T] but for the B* = 400 Gauss dipole 
case. The maximum velocity vector corresponds to 190 km s~^. 

estimate of Paper I and accounts for approximately 30% 
of their overestimate of ta- 

As described in M att fc Balicld ()2004l ). a given sim- 
ulation can be scaled to other systems with the same 
characteristic velocity ratios, so that the resulting value 
of taZ-R* is valid for a family of solutions. The sim- 
ulation with / = 0.004 scales to a solution very sim- 
ilar to the solar wind with i?* — IRq, A/, — IMq, 
Afw = 1-3 X IQ-'^^Mq yr-\ = 1.5 G, and all speeds 
are increased by a factor of 2. Thus, for these parame- 
ters, this simulation predicts a lever arm length of 8.33 
Rq for the case of the solar wind, and t„ — 6.8 x 10^^ 
er g. This torque is consisten t with the numerical results 
of lWashimi fc Sakurail (|1993l ). but a factor of a fe w time s 
smaller than observationally determined values ()Lilll999f ). 
To obtain the observed solar torque, corresponding to 
rA = 12.2i?Q, the simulation would require (e.g.) a sub- 
stantially str onger mag netic field than 1 G. This was also 
suggested by |lJ (|1999f ) and our simulations corroborate 
that suggestion. As the estimate in Paper I assumed the 
canonical value of _B, w 1 G for the sun, this accounts for 
most of the discrepancy between our simulation results 
and the Paper I estimate of ta- 

To capture a range of spins appropriate for the T Tauri 
star "slow rotators," we also ran cases with spin rates of 
twice and half of the fiducial spin rate. The results of 
these simulations are listed in the 3rd and 4th row of 
table H 

4.2.2. Dipole Field Strength 

Measurements for th e mean \Bf,\ exist f or a number 
of T Tauri stars (e.g., Johns- Krulll l2007b[) . These re- 
sults show a remarkably consistent field strength for all 
stars of around 2 kG. Measurements of the longitudi- 
nal field (which limits the global, dipole c omponent) ex- 
ist only for a handful of accreting stars (iBouvier et al.l 
[2007t l,Tohns-K7^[2n07ah . These measurements suggest 
the dipole component is no greater than 200 G (though 
larger values would be allowed for special viewing ge- 



FlG. 4. — Same as figure [T] but for the 2 kG quadrupole case. 
The maximum velocity vector corresponds to 170 km s~^. 

ometries). Given the small number of measurements, it 
is still relevant to consider stronger dipole field strengths. 
Thus we have run cases with i?, — 400 and 2 kG. 

Figure [3] illustrates the wind solution for the case with 
i?* = 400 G, and the results of both cases are listed in 
table [2] It is clear that the strength of the field has a 
strong influence on the stellar wind torque. 

4.2.3. Surface Field Geometry 

The fact that T Tauri stars have a mean field of 
|-B*| ~ 2 kG with a much weaker dipole component, indi- 
cates that the stellar surface field is dominated by higher 
order multipole fields. Therefore, it may be important 
for future work to include much more structured fields 
than we consider here. To begin to quantify the effects of 
higher order fields, we ran two cases that were initialized 
with a quadrupolar field of the form 

BR = B,{?,cos'e~l) {?^^ 

Be = 2B^ cos 9 sin 9 (12) 

in spherical coordinates. 

We ran cases with equal to 1 kG and 2 kG, listed 
as "1 kG quad." and "2 kG quad." respectively in ta- 
ble [21 Figure [5] illustrates the wind solution for the 2 
kG quadrupole case. It is clear from the figure that, 
compared to the cases with a dipole field, the shape of 
Alfven surface is quite different. Also, for a given value 
of B^,, the effective lever arm length is much shorter for 
the quadrupole case. 

It is evident from table [2] that the stellar wind torque 
from a star with a 200 G dipole field is comparable to that 
with a 1-2 kG quadrupole field. Thus, the stellar wind 
torque is very sensitive to surface field geometry. How- 
ever, since the measured surface field strengths of ^ 2 kG 
likely include contributions from even higher multipoles 
than a quadrupole, it seems that the dipole component 
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will generally dominate near the Alfven surface. So the 
strength of the dipole component should generally be the 
most important for determining the torque. 

4.2.4. Mass Loss Rate 

We have thus far considered quite massive winds, moti- 
vated by recent suggestions in the literature for accretion 
powered winds. However, the value of is very uncer- 
tain and is likely to exhibit a wide range in values from 
one object to the next. In addition, it would be inter- 
esting to predict what torques may be expected for the 
winds from the non-accreting, weak line TTSs. We ex- 
pect these stars to have solar-like winds that are quite 
enhanced relative to their main sequence counterparts, 
yet probably less powerful than winds from the accreting 
stars. 

Unfortunately, our method is limited to cases with 
lever arms that are not too long, since longer lever arms 
requires a larger Alven speed on the stellar surface. A 
large Alfven speed increases the time for the simulation 
to run and also increases the error in the solution (e.g., 
by increasing the effective diffusion rate in our code) . For 
this practical reason, we were limited to running only two 
cases with lower covering a range in of a factor 
of 50. These are listed in the 9th and 10th rows of table 

m 

4.2.5. Stellar Radius 

T Tauri stars contract as they age, so stars of a given 
mass exhibit a range of radii during this phase. Thus, 
it is important to consider here different combinations of 
stellar mass and radius. Table [2] contains results from 
two cases with i?, = l.Si?© and i?* = 3Rq. Note that 
changing changes Uesc, so these cases have a different 
coronal temperature and fJ*, in order that Cg/wosc and 
/ are constant. From the values in the table, it is evi- 
dent that the stellar wind torque is very sensitive to i?, . 
The reason for this is twofold. First, since i3* is fixed, 
a larger stellar radius corresponds to a larger dipole mo- 
ment (/X = B^i?^), which is capable of conveying a larger 
torque. Second, a larger stellar radius decreases the sur- 
face gravity, and so the influence of the magnetic field 
relative to gravity is increased (i.e., wa/wcsc increases). 
Thus, {rxlR*Y increases with i?,, and though fi, de- 
creases (to keep / fixed), the quantity increases, 
so the net torque increases. 

4.2.6. Stellar Mass 

Cases with half and twice the fiducial stellar mass are 
also listed in table [21 As with the cases of different i?, , 
note that changes in change Wcsc, so we have adjusted 
the coronal temperature and fi, to keep the parameters 
listed in table [1] fixed. As with the case of varying i?,, 
a change in Vesc changes the relative importance of the 
magnetic field with the gravity. Thus, va/R* is larger 
for a smaller M*. However, since we have fixed /, a 
smaller M, means a smaller f2» so that the net stellar 
wind torque decreases. 

4.2.7. Wind Acceleration 

The increase of the wind speed with distance from 
the star depends on the details of the wind acceleration 



mechanism. In a Parker wind, the temperature (param- 
eterized by Cs/wosc) and the cooling/heating of the gas 
as it flows (parameterized by 7) are the key physical 
properties determining the velocity profile in the wind. 
A hotter wind accelerates more rapidly and achieves a 
higher speed than a cooler wind. A wind with a larger 7 
(closer to 5/3) cools more rapidly as it expands, and so 
the bulk of the acceleration takes place closer to the star. 
Similarly, if the wind is instead accelerated by something 
other than thermal pressure, the velocity profile may be 
altered. 

In order to quantify the effect of varying the accelera- 
tion in the wind, within the framework of the pressure- 
driving mechanism used here, we have run three more 
simulations. The results of these are listed in the last 3 
rows of table [21 In order, these represent winds that are 
hotter, colder, or with less heating (i.e., more adiabatic 
cooling) than the fiducial case. The relatively large effect 
of Cs/fesc and 7 on the wind speed near the stellar sur- 
face is evident by the very different values of required 
to keep fixed, listed in table [H The wind velocity 
at the base of the corona, for fixed M^j, varies as the 
inverse of the variation in So these cases represent 
large differences in the wind acceleration rate. The ef- 
fect on the torque is relatively small, but is not entirely 
negligible. It will be important for future work to deter- 
mine the wind torques for different driving mechanisms. 
The preliminary conclusion to be drawn from this work is 
that the wind velocity profile, and therefore wind driving 
mechanism, does not have a large effect on the torque. 

5. SEMI- ANALYTIC FIT FOR THE EFFECTIVE ALFVEN 
RADIUS 



In section [2711 we pointed out that no reliable formula- 
tion exists for predicting the Alfven radius (and therefore 
torque) in a stellar wind from fundamental parameters. 
However our parameter study, even though somewhat 
limited, can be used to provide a numerically based ap- 
proach to this question. We will use the result future 
work, and it also will be of general interest for other stel- 
lar wind studies. 

In a 1-dimensional theory, one can assume that the 
magnetic field strength approximately follows a single 
power law of the form B — B^,{R^./ R)" . Then the condi- 
tion that t he wind speed equals the Alfven spe ed at ta 
gives fe.g.. lKawale"i1ll988HTout fc Pringlelll992f ) 



R* 



BiR^ 



(13) 



rA 



where VrA is the wind speed at the Alfven radius. There 
are a number of problems. First, the true magnetic field 
stren gth in a wind does no t follow a single power law 
(e.g.. lMestel fc Spruij|1987[) . Second, the Alfven surface 
is neither a sphere nor a cylinder and a spherical model 
is quite misleading. Third, and perhaps most vexing, is 
the fact that VrA has different values at different points 
along the Alfven surface and cannot be determined a 
priori. Finally, there is no explicit dependence of ta on 
the spin rate or driving properties of the wind, which we 
also know to be false. 

There is, however, a more general way of scaling the 
Alfven radius that is suggested by basic theory. Another 
clue is that since the numerical simulations are carried 
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Fig. 5. — Effective lever arm le ngth in the stellar wind versus the 
quantity in brackets in equation II14I I. Shown are the results of our 
entire parameter study including the fiducial case (filled circle); 
cases with different spin rates (pluses); cases with a quadrupole 
field (triangles); cases with different Cb/v^bc or 7 (squares); and 
all other cases (diamonds) , representing those with different values 
of Bt,, Rt, Mvi, or M, . The line represents the best fit to the 
fiducial and "other" cases, given by equation I I14I I with K 2.11 
and m 0.223 

out in normalized units, they are scalable to any system 
with with similar characteristic speeds on the stellar sur- 
face fe.g.. iMatt fc Balick|[200l . This suggests that we 
can replace VrA with the stellar surface escape speed and 
calculate the Alfven radius using 



R* VMwWesc 



(14) 



where K and m are dimensionless constants. The 
quantity inside the bracket measures the effective mag- 
netization of the wind, is similar to that used by 
lud-Doula fc Owockil (|2002l see their eq. 7), and arises 
naturally in disk wind theory (e.g., equation 2.27 of 
iPelletier fc Pudrit3ll992l ). It is also a quantity that can 
be fixed by observations of stellar properties and wind 
parameters. 

In figure [5l we plot ta/R* as a function of the quantity 
in brackets in equation (fT4|) on a log-log scale, for all 17 
of our simulations. Since we know equation (jl4p does not 
properly include the effects of stellar rotation or the wind 
driving mechanism, we calculate the best fit K and m to 
the fiducial case and only those cases with variations on 
B.,, i?*, Mw, and M,. The fit, giving K w 2.11 and 
m « 0.223, is plotted as a line in the figure. 

It is remarkable that this fit matches all of the relevant 
simulations (filled circle and diamonds in the figure) to an 
accuracy of less than one percent. This is at the level o f 
precision of the numerical method (|Matt fc Bahckl[200^ . 
Remember that we have taken Cg/wosc and / as our pa- 
rameters, so that cases with a different value of Vcsc (i-e., 
those with different i?* and Af») actually also have differ- 
ent wind temperatures (i.e., Cg) and stellar angular spin 
rates (f2*) than the fiducial case. If we had chosen il* 



and Cg as our fixed parameters, there would be a lot more 
scatter of the diamonds around the line in figure [5l Fur- 
thermore, our simulations self-consistently capture the 
interaction between the stellar wind, magnetic field, and 
rotation, without resorting to assumptions about (e.g.) 
the extent of the dead zone, the dependence of magnetic 
field strength with radius, or latitudinal variations in 
wind quantities. Therefore, our semi-analytic formula- 
tion appears to be an improvement over existing theory. 

Equation ITU) does have some limitations. Neither the 
previous analytic formulations nor our own semi-analytic 
approach properly includes the effects of varying stellar 
rotation or wind driving (as evident in figure [5]). As an 
illustrative example, in a smaller parameter study with a 
rotation rate comparable to the solar rate (not presented 
here), we found that X « 3.0 and m « 0.19. Also, note 
that a line connecting the two points with a quadrupole 
field suggests K « 1.7 and m « 0.15, for these cases. 
Thus, the indicies K and m are quite sensitive to the 
field geometry and have a smaller (but non-negligible) 
sensitivity to changes in the stellar spin rate and the 
wind acceleration rate/mechanism. We leave the precise 
determination of the sensitivity of ta to these parameters 
for future work. 

We can now combine equations ([3]) and ([M]) to get a 
formula for the stellar wind torque. 



^ f 1- 

V2 ^ """^ 



2m 



M. 



1 — 2m ^l-|-4m 



B 



4m 



(15) 



though we know this does not properly contain the de- 
pendence (e.g.) on /. This equation is essentially the 
same as that derived bv iKawaleil ()1988l ). except for the 
value of the the dimensionless constant out front and of 
the expected value of the exponent parameter m. The 
constant is not so crucial, and usually this can be cali- 
brated to the solar wind torque for a predictive theory 
(though we have not done this here). On the other hand, 
the value of m is of far greater importance for predicting 
the torque for a range of parameters. 

In particular, authors typically have chosen a power 
law such that the stellar wind torque is nearly or 
completely independent of Af^ (effectively, 
0.5), which results in r, 



v..^,, .^^^.^^ ... ,w (X i?f (e.g.. iKawalej 

19881 iPinsonneault et all 119891 : iBarnes fc Sofial Il996t 



Bouvier et al.lll997D . Ou r basic u nderstanding of the ob- 



served Skumanich-style (jSkuman ich 1972) spin down of 
main sequence stars (fi* cx <~^/^), as well as the expected 
dependence of magnetic field strength with rotation rate 
(jBelcher fc MacGregoij [l976l ). appears to rely on this 
or a similar formulation. Our fit value of m w 0.223 
gives approximately Tw o c i?^'^ (this was also found by 
IWashimi fc Shibatalll993f l. which is substantially differ- 
ent. 

In order to understand the difference between our value 
of the exponent and that used by others, it is instructive 
to consider the power law formulation of the magnetic 
field used to derive equation (fT3l) . It is generally expected 
(e.g.. lMestelll98l iMesFel fc SDruitlll987l : lKawalei)[T98ah 
that when the surface magnetic field is dipolar, as we 
are considering here, the effective power law index of the 
magnetic field in the flow will lie somewhere between the 
value for a dipole (n = 3) and that for a split monopole 
(n — 2). This has been the primary justification for the 
power laws used in the literature. Indeed our simulations 
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display the expected behavior of exhibiting a dipolar ge- 
ometry near the star and an approximately monopolar 
geometry far from the star (e.g., figure [1]). However, by 
comparing equations ([13]) and (fT4| . we see that our fit 
value of m « 0.223 seems to imply a magnetic power law 
of n« 3.2. 

It is important to realize that the divergence of the 
magnetic field in the flow, captured by the power law 
index n, is not the only important effect, and this is 
why the formulation of equation (|13p is misleading. Here 
are two reasons. First, using 1-dimcnsional reasoning, 
in an accelerating wind, the behavior of u^a mitigates 
the response of to the parameters. For example, for 
an increase in B^, the Alfven radius will become larger, 
but since the flow is accelerating, VrA will also increase. 
iKawaleil ([l988, ) made the approximation that VrA equals 
the escape speed at ta. In this case, VrA decreases with 
radius, giving the opposite effec t of an a ccelerating wind. 
Similarly, the approximations of lMestel (jT984) that u^a is 
constant for a slow rotator and proportional to Q^r^ for 
a fast rotator do not well-approximate the acceleration 
exhibited in the winds we simulated. The second reason 
for the surprisingly weak dependence of ta on parameters 
is in the amount of open magnetic flux that participates 
in the flow, which again is not included in the derivation 
of equation (|13p , and which again mitigates the effect of 
parameters on ta- For example, for an increase in B^, 
a smaller area on the stellar surface wi ll have open flux 
(e.g., compare figures [1] and [3l and see iMestel fc Spruit! 
Il987( ) , so TA will not increase as much as expected in the 
magnetic power law formulation. 

In future work, it will be important to extend equation 
([14]) to include the effects of rotation, etc. Furthermore, 
much work is needed to explore the full consequences 
(e.g., for main sequence stars) of the significantly smaller 
exponent we find, compared to many previous works. 

6. SUMMARY AND CONCLUSIONS 

Using 2D (axisymmetric) MHD simulations, we com- 
puted steady-state, stellar wind solutions for a parame- 
ter range appropriate for T Tauri stars. We carried out a 
parameter study including variations of the stellar mass, 
radius, surface magnetic field strength, and rotation rate, 
as well as mass loss rate, wind acceleration rate, and two 
different magnetic geometries (dipole and quadrupole). 
Our solutions enabled us to determine the angular mo- 
mentum carried in the wind, and its dependence on many 
of the parameters of the system. Our main conclusions 
can be summarized as follows: 

1. For fiducial parameters, the torque is of the same 
order (~ 10'^^ erg) as estimated in Paper I. There- 
fore, if the stellar winds of TTSs have similar pa- 
rameters to those considered here, they should have 
a significant influence on the stellar spin. 

2. The stellar winds are in the regime of moderately 
fast magnetic rotator winds. They produce jets, as 
well as a wide-angle flow (see, e.g.. lMatt fc Bailed 
l2004f ). which should interact with, and be modi- 
fled by, surrounding material (not included in our 



simu lations; e.g.. lGardiner et al]|2003l : IShang et al.l 
l200fil ). 

3. The cases with quadrupole fields resulted in a 
torque that is much weaker than cases with a dipole 
field of the same surface field strength. Specifi- 
cally, we find that a 200 G dipole field exerts the 
same stellar wind torque upon a star as a 1-2 kG 
quadrupole. This illustrates the very strong effect 
of magnetic geometry on the stellar wind torque. 

4. We ran cases where the mass loss rate and other 
parameters were fixed, but the thermal wind driv- 
ing parameters were varied. For large variations in 
the wind acceleration, the torque changed by less 
than a factor of 2. This suggests that the details 
of the velocity profile are not of fundamental im- 
portance, and our solutions should be a reasonable 
approximation for winds with other wind driving 
mechanisms. However, it will still be important 
for future work to compare our torque results to 
stellar wind solutions that use alternative driving 
mechanisms. 

5. Our determination of the torque allowed us to cal- 
culate the Alfven radius (via eq. [3]) , which is a fun- 
damental quantity in MHD wind theory. We com- 
pared our numerical solutions to previous analytic 
work and obtained a semi-analytic formulation for 

rA/R* (X [Bji?^/(MwWosc)]'", with m w 0.22 (eq. 
[T4| . that well-describes many of our simulations 
with dipole fields. This formulation appears to be 
an improvement over existing work, and the ex- 
ponent m is significantly smaller than usually as- 
sumed. 

We will continue to develop the theory of accretion- 
powered stellar winds in forthcoming work. In a compan- 
ion paper (the third in our series), we compare the stellar 
wind torques computed here to the torques expected to 
arise from the interaction between the star and an ac- 
cretion disk. We find spin-equilibrium (net zero torque) 
solutions and test the suggestion of Paper I. In a later 
paper, we will use the stellar wind solutions of this work 
to compute emission properties of TTS coronal winds. 

Initial results of various aspects of this work were pre- 
sented at a number of meetings, and we wish to thank 
many participants who provided interesting questions, 
discussion, and ideas including: Gibor Basri, Sylvie 
Cabrit, Steve Cranmer, Andrea Dupree, Suzan Edwards, 
Christian Fendt, Will Fischer, Shu-ichiro Inutsuka, Chris 
Johns-KruU, Marina Romanova, Frank Shu, Keivan Stas- 
sun, Asif ud-Doula, Jeff Valenti, and others. We also 
thank the referee, Ruben Krasnopolsky, for his useful 
suggestions for improving the paper. SM is supported by 
the University of Virginia through a Levinson/VITA Fel- 
lowship partially funded by The Frank Lcvinson Family 
Foundation through the Peninsula Community Founda- 
tion. REP is supported by a grant from NSERC. 
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